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摘 要 ”本文 将 收敛 速度 快 的 变形 玻 恩 迭代 法 (DBIM) 结合 数值 模式 匹配 法 (NMM) 用 于 实际 
工程 中 的 阵列 感应 成 象 仪 (AIT) 的 测量 信号 的 反 演 和 剖面 成 象 . 在 反 演 迭代 过 程 中 , 利用 NMM 
解 的 z 向 解析 性 质 ， 避 免 了 双重 数值 积分 ， 提 高 了 计算 效率 ， 同 时 也 提高 了 反 演 质量 。 文 中 利 
用 不 完备 的 测量 数据 对 AIT 进行 了 反 演 ， 表 明 该 方法 非常 适合 实际 工程 应 用 . 

关键 词 ” 逆 散射 ， 阵 列 感应 成 象 仪 ， 变 形 玻 恩 欠 代 ， 数 值 模式 匹配 法 
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1 引言 


电磁 逆 散 射 是 当今 电磁 学 研究 的 热点 之 一 . 它 已 被 广泛 应 用 于 地 球 物理 勘探 、 医 学 成 
象 、 目 标识 别 以 及 材料 无 损伤 检测 等 实际 工程 领域 中 . 仅 以 地 球 物理 应 用 而 言 ， 近 年 就 吸引 
了 众多 学 者 的 研究 上 L- 有 .众所周知 ， 收 、 发 线圈 系 之 间 有 很 强 的 直 耦 信号 ， 要 在 这 么 强 的 
育 景 场 中 提取 散射 场 信息 是 非常 困难 的 ， 所 以 在 实际 工程 中 为 了 对 消 直 耦 信号 通常 采用 复 
合 线圈 系 . 实际 工程 中 常用 的 阵列 感应 成 象 仪 (AIT) 就 是 根据 几何 因子 理论 , 用 软件 聚焦 的 
方式 来 提高 其 纵 问 分 辨 率 和 增加 径 向 探测 深度 . 由 于 几何 因子 理论 是 建立 在 均匀 介质 模型 
基础 上 ， 又 没 考虑 非 线性 部 分 ， 所 以 是 一 种 近似 方法 ， 其 解释 效果 必然 受到 了 限制 . 

本 文 就 是 针对 AIT ， 基 于 严格 的 积分 方程 出 发 ， 在 非 均 匀 背 景 介 质 中 采用 了 变形 玻 恩 
达 代 法 (DBIM) 来 重建 二 维 轴 对 称 介 质 结构 的 电导 率 分 布 。 由 于 该 方法 在 反 演 迭 代 过 程 中 
需要 多 次 调用 正 演算 法 ， 所 以 选用 一 种 高 效 精确 的 正 演 方法 是 提高 反 演 质量 和 计算 效率 的 
关键 ,本 文中 正 演 方法 选用 了 高 效 的 数值 模式 匹配 法 NMM) . 由 于 NMM 解 在 z 向 的 解析 
性 质 ， 故 其 在 反 演 过 程 中 z 向 积分 有 解析 解 ， 加 之 径 向 用 线性 插值 方法 ， 和 避免 了 双重 数值 积 
分 ， 与 采用 数值 积分 方法 罗 相 比 ， 其 计算 效率 提高 了 一 个 数量 级 . 


2 非 线性 积分 方程 的 建立 和 变形 玻 恩 迭代 法 


企 二 维 轴 对 称 非 均匀 介质 中 ， 假 设 磁 导 率 J(p,z) = ko ， 则 只 有 电导 率 o = o(p,z) 和 介 
电 贡 数 e= e(p,z) 是 pp 和 > 向 函数 . 在 上 述 条 件 下 ， 位 于 轴线 上 细 环 天 线 仅 产生 Es 分 其 的 
场 ， 其 散射 场 满足 以 下 积分 方程 是 : 


sca 1 "I ! ne / minc 1 .1 1 .1 1 1 
Es (pz pt, zt) iwnuol 局 dz / dp Eo (p,zZ,p ;之 )0O(O ;之 )Eo(p ;2 , Pt, 2t). (1) 
上 式 右 端 Bo(p',> pb 和 ) 也 是 目标 函数 60(p',z') 的 函数 ， 所 以 (1) 式 是 非 线性 积分 方程 . 
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Rs R, Rs R; RT R, R, Re 
图 1 阵列 感应 成 象 仪 示意 图 


由 于 我 们 所 考虑 的 是 实际 工程 中 所 应 用 的 复杂 线圈 系 一 -阵列 感应 成 象 仪 (AIT) ， 它 是 
由 八 对 三 线圈 系 基 本 单元 组 成 的 ， 每 单元 是 由 一 个 发 射 天 线 和 两 个 绕 向 相反 经 直 耦 信号 对 
消 的 接收 天 线 构成 ， 如 图 1 所 示 . 

对 于 AIT 的 三 线圈 系 单 元 ， 各 接收 线圈 对 上 的 接收 场 应 为 


Em (p2, 22, p1, 21, pt, Zt ) 一 No1 Es (po2, 22, pt, zt) Egy(p1, 21 pt, Zt)) (2) 


式 中 N21 为 线圈 对 的 臣 数 比 。 由 (1) 式 和 互 易 定 理 可 将 上 式 写 为 
Em(p2, 2, pl1, 1, pt, Zt) [N21 ES (p2, 22, pt, zt ) 2 Ey°(p1, z1, pt, 2t)] 


十 CO 十 oo 
= |/ dz |/ dp'[N21 Eg (p',z ,Pp2, 22) bs Eg(p',z ,pi1,21)]60(p',z ) Eo(p', 2 , pt, zt)- (3) 
一 Co 0 


在 弱 散 射 情形 下 可 应 用 玻 恩 近似 ， 令 积分 号 中 的 Bo (pz pi, 和) 守 Be(p',z', pb zt) ， 
而 入 射 场 只 与 背景 介质 有 关 ， 与 50(p',z') 无 关 ， 故 将 (3) 式 积分 方程 线性 化 . 
为 求 出 未 知 目标 区 的 电导 率 , 需 把 目标 区 的 剖面 分 布 离散 化 , 并 按 脉冲 基 郴 数 序列 {f;(p,z)} 
离散 (3) 式 ， 得 到 如 下 和 矩阵 方程 : 
e=M.a (4) 
式 中 e 表示 测量 数据 矢量 (Nn 个 元 素 ) ，a 是 指 目标 函数 矢量 (Nm 个 元 素 ) ，M 为 Nn x Nm 
的 矩阵 ， 其 官 阵 元 素 表示 为 


十 co 十 Co 
Min,; = a az | dp’'f;(p',z' )[No1 Eg(p', 2', p2l, Z21) 
一 Es(p',， Ss pll) zu] Es (p', 2 Dv Sin) (5) 
上 式 中 mn 和/ 分别 表示 第 n 个 发 射线 圈 和 第 /个 接收 线圈 对 . 
反 演 问题 也 就 归结 为 求解 (4) 式 中 的 a ， 由 于 (4) 式 具 有 不 适 定性 ， 通 过 Tihonov 正则 
化 方法 将 其 变 为 适 定 问题 。 首 先 应 构造 一 代价 函数 : 
Cla)= le— Mall +)R .oll, (6) 
式 中 五 为 正则 化 矩阵 (文中 应 用 单位 矩阵 ) ， 7 为 正则 化 参数 . 要 使 (4) 式 有 稳定 的 解 ， 那 
么 (6) 式 的 范 数 应 最 小 ， 即 变 分 为 零 ， 由 此 可 得 
a= (VM' .M+7yR .RBR-!.M .e. (7) 
式 中 + 表示 复 共 轿 转 置 . 本 文采 用 共 辆 梯度 法 求解 该 方程 便 可 求 得 a . 在 第 一 次 DBIM 中 ， 
首先 对 a 赋 一 初 值 ， 通 过 解 (7) 式 中 的 a 去 更 新 背景 介质 中 的 电导 率 分 布 ， 在 新 的 背景 介 
质 中 计算 背景 介质 的 入 射 场 B82° ， 再 通过 解 (7) 式 又 可 得 到 新 的 电导 率 分 布 ， 迭 代 过 程 进 
行 到 散射 场 的 计算 值 和 测量 值 之 间 的 相对 剩余 误差 (RRE) 小 于 给 定 值 时 迭代 结束 ， 此 时 求 
得 的 电导 率 分 布 便 是 介质 结构 的 电 参 数 分 布 的 重建 结 
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3 高效 的 正 演 方法 一 一 数值 模式 匹配 法 


利用 DBIM 进行 反 演 时 ， 每 次 迭代 均 要 调用 正 演 程序 . 正 演 计算 本 文 甩 用 了 高 效 的 
NMM ， 各 区 的 场 由 文献 [6] 可 知 ， 其 中 源 区 (p) 的 场 表示 为 


ci 1 之 一 1; zz 一 1 f (p") i 
-fp(P)[e 和 + @ plr dp)G, pr1e 有 (之 之 为 
peg wl | ， f (p") ® 
-Oe + et Gye te Boi De (2 <) 
在 源 上 方 区 域 的 场 (n > p) 中 ， 
i 
pBne = — F(t + either-t) Ginter tt) An det. (09) 
2 ks,: Do» 


在 源 下 方 区 域 的 场 (m < p) 中 ， 


f,(p') 
ks Ds 


pBmo = = fh (Pleam 7s) + tems (edn) Gn nemr (dn)] Br, (10) 
上 式 中 4 和 B 为 场 幅 ， G 为 广义 反射 系数 矩阵 ， 它 计 入 界面 的 多 次 反射 效应 f(p) 则 
表示 第 n 层 的 径 向 本 征 模 ， 径 向 格 点 间 的 本 征 值 由 线性 插值 得 到 ， 即 

fn(p) = krp + br. (11) 


4 二 重 积分 的 解析 解 


由 (5) 式 可 知 ， MM 矩阵 元 素 是 由 两 个 场 乘 积 的 二 重 积分 得 到 .我 们 知道 计算 二 重 数值 
积分 是 非常 繁 锁 的 ， 尤 其 是 在 积分 区 域 比较 大 以 及 计算 精度 要 求 较 高 时 ， 使 计算 量 显著 增 
加 ， 甚 至 计算 机 都 无 法 运算 . 径 向 的 数值 本 征 模 经 线性 插值 后 也 可 表示 成 直线 方程 ， 其 径 向 
积分 非常 容易 。 为 讨论 方便 我 们 假定 收 、 发 线圈 均 在 目标 区 下 方 ， 场 由 (9) 式 来 表达 ， 那 么 
z 问 积 分 可 写 为 下 式 : 


dn 
| Eg(p',z', pri, Zrl) Bg(p',z', pin, ztn)dz = M; + Mi + M+ MS, (12) 
dn—i1 
其 中 
Mi = 元 二 一 二 [eiti(dn dn -1) eikr(dn dn-1) — 1.0]; (13) 
t rT 
三 | 
ii 十 Br)(dn 一 dm -1) _ ,2ikn(dn—dn-1) 5 
[到 € 9 天 fi 
a | fy (14) 
(dn 一 UR ks = kr.; 
1 2iki(dn—dn-1) i(ktt+kr )(dn—dn-1) 。 
— le 总 沁 k Kk 
me (15) 
(ed et a ks = kr; 
MS = ? [本 人 Sa (16 ) 


ki 十 kr 


> 


电导 率 (S/m) 
be 
电导 率 (S/m) 


(b) 
0.020 门 
by 
0.015 
8 仿 
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图 2 不 完备 测量 信息 的 径 向 反 演 结果 : (a) 原 电导 率 分 布 ; (b) 第 一 次 DBIM 友人 代 绪 采 
(c) 第 八 次 DBIM 达 代 结果 ; 〈d) 场 的 相对 剩余 误 老 


太 和. 是 目标 区 的 本 征 值 其它 几 种 情况 也 可 类 似 得 到 .此 时 在 计算 M 矩阵 时 就 避免 了 
二 重 数值 积分 ， 大 大 提高 了 计算 效率 ， 在 大 区 域 反 演 时 效果 尤为 明显 ， 


5 数值 计算 结 采 及 结论 


在 以 下 例子 中 , 背景 介质 均 取 为 非 均 匀 介 质 结 构 ， 即 井 外 背景 介质 电导 率 和 并 内 电导 率 
分 别 为 0.01S/m 、 0.1S/m ， 所 使 用 频率 为 20kHz ， 对 图 1 所 示 的 阵列 感应 成 象 仪 的 测 基 
数据 进行 反 演 .由 于 仪器 在 井下 是 沿 井 轴 移 动 测 基 ， 对 径 向 剖面 而 言 其 测 基数 据 是 不 完备 
的 ， 实 际 工 程 应 用 中 通常 是 串 接 径 向 探测 深度 不 同 的 仪器 组 合 使 用 , 希望 获得 不 同 径 向 深度 
的 电导 率 , 即使 是 这 样 也 很 难 获得 满意 的 结果 . 现在 我 们 来 考查 一 下 仪器 的 径 向 和 纵 回 分 辨 
率 . 图 2(a) 是 两 个 0.25mx0.25m 相距 0.25m 的 真实 电导 率 分 布 , 图 2(b) 是 DBIM 第 一 次 迭 
代 结 果 ， 它 还 分 辨 不 出 第 二 个 介质 块 ， 图 2(c) 是 DBIM 第 八 次 迭代 结果 ， 它 已 经 非常 清晰 
地 分 辨 出 两 个 介质 块 . 图 3(a) 是 两 个 0.15mx0.75m 相距 0.15m 的 真实 电导 率 分 布 ， 图 3(b) 
是 DBIM 第 一 次 迭代 结果 ， 其 分 辨 的 效果 并 不 好 ， 图 3(c) 是 DBIM 第 八 次 欠 代 绪 果 ， 它 已 
经 能 够 很 好 地 分 辨 出 两 个 介质 块 ， 而 且 其 最 大 值 已 接近 真 值 。 上 面 两 个 例子 说 明 ， 用 DBIM 
方法 对 AIT 成 得 具有 很 高 的 径 向 和 纵向 分 辨 率 ， 而 用 传统 的 软件 聚焦 的 方法 是 达 不 到 这 样 
的 效果 . 

本 文 将 变形 玻 因 过 代 方法 并 结合 数值 模式 匹配 理论 应 用 于 阵列 感应 测 井 的 成 象 。 由 于 
NMM 是 半 解 析 半 数值 的 方法 ， 利 用 其 解 的 z 向 解析 性 质 可 获得 z 向 积分 的 解析 结果 ， 而 在 
径 向 格 点 间 我 们 采用 了 线性 插值 的 方法 ， 从 而 避免 了 积分 方程 中 的 二 重 数值 积分 . 与 数值 积 


1 期 


杨峰 等 ， 基 于 积分 方程 的 阵列 信号 的 反 演 方法 


上 一 
一 


电导 率 (S/m) 
© 


un 


1 > > 


电导 率 (S/m) 


se 


一 一 
= 


DD 
Dem 


二 让 


| 
一 
谋 me 


之 RS 
(a (b) 
0.025 
灿 0.020 
下 张 
代 0.015 
方 一 
& 0 区 0.010 
由 到 
这 四 0.005 
> 6 7 
迭代 次 数 
ES (d) 


不 完备 测量 信息 的 纵向 反 演 结果 : ”(a) 原 电导 率 分 布 ; (b) 第 一 次 DBIM 选 代 结果 ; 


(c) 第 八 次 DBIM 和 迭代 结果 ; (d) 场 的 相对 剩余 误差 
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分 相 比 ,这 样 做 的 好 处 是 计算 量 显 著 减少 ,特别 是 对 于 大 区 块 齐 面 的 反 演 更 显示 出 其 明显 的 
优越 性 ， 如 利用 Pentium 133 微机 在 NDP-FORTRAN V3.20 环境 下 ， 本 征 函 数 取 52 个 ， 对 


图 2 所 示 5x5 离散 块 ， 56 个 测试 数据 ， 迭 代 一 次 计算 时 间 不 超过 10min ， 


对 厚 层 和 大 区 域 


进行 反 演 成 象 ， 计 算 工 作 量 并 不 增加 .通过 上 面 的 例子 表明 ， 经 过 适当 的 和 迭代 次 数 后 ， 反 演 
结果 已 十 分 接近 真实 电导 率 分 布 。 此 方法 非常 适合 实际 的 测 井 解释 工作 应 用 . 
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INVERSION OF THE MEASURING DATA OF 
ARRAY INDUCTION IMAGER TOOL BASED ON 
THE INTEGRAL EQUATION 


Yang Feng Nie Zaiping Zhao Yanwen Zhang Yerong”* 


(Dept. of Microwave-Engin., Univ. of Electro. Science and Tech. of China, Chengdu 610054) 
“(Dept. of Electro. Engineeriny, Huazhong University of Science and Technology, Wuhan 430074) 


Abstract This paper deals with the profile reconstruction of the Array Induction Imager 
Tool( AIT) measurement using the Distorted Born Iterative Method(DBIM) and an efficient Nu- 
merical Mode-Matching( NMM) method for the conductivity distribution in a two-dimensional 
axisymmetric medium. In each iteration of DBIM, the NMM method is applied to solve a 
forward solver. Meanwhile, the analytic expression of the integration in the nonlinear integral 
equation has been derived and then the inversion of the efficiency and accuracy is improved dra- 
matically. Several examples have shown that this method is practicable for the interpretation 
of electrical logging data. 


Key words Inverse scattering, Array induction imager tool, Distorted Born iterative method, 
Numerical mode-matching method 
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